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ABSTRACT 

Cold Dark Matter (CDM) has become the standard modern theory of cosmo¬ 
logical structure formation. Its predictions appear to be in good agreement 
with data on large scales, and it naturally accounts for many properties of 
galaxies. But despite its many successes, there has been concern about CDM 
on small scales because of the possible contradiction between the linearly rising 
rotation curves observed in some dark-matter-dominated galaxies vs. the 1/r 
density cusps at the centers of simulated CDM halos. Other CDM issues on 
small scales include the very large number of small satellite halos in simulations, 
far more than the number of small galaxies observed locally, and problems con¬ 
cerning the angular momentum of the baryons in dark matter halos. The latest 
data and simulations have lessened, although not entirely resolved, these con¬ 
cerns. Meanwhile, the main alternatives to CDM that have been considered 
to solve these problems, self-interacting dark matter (SIDM) and warm dark 
matter (WDM), have been found to have serious drawbacks]]] 


1 This paper is a significantly updated revision of 



1 Introduction 


Sometimes a theory is proposed in relatively early stages of the development 
of a scientific held, and this theory turns out to be not only a useful paradigm 
for the further development of the held it also survives confrontation with 
a vast amount of data, and becomes accepted as the standard theory. This 
happened with General Relativity , and it seems to be happening now with 
general relativistic cosmology. It appears that the universe on the largest scales 
can indeed be described by three numbers: 

• Hq = 100/ikms -1 Mpc -1 , the Hubble parameter (expansion rate of the 
universe) at the present epoch, 

• il rn = p/p c , the density of matter p in units of critical density p c = 
3-ffg(87 tG) -1 = 2.78 x 10 n h 2 M e Mpc” 3 , and 

• Q,\ = A(3I7q) - 1 , the corresponding quantity for the cosmological con¬ 
stant. 

The currently measured values of these and other key parameters are summa¬ 
rized in the Table below. It remains to be seen whether the “dark energy” rep¬ 
resented by the cosmological constant A is really constant, or is perhaps instead 
a consequence of the dynamics of some fundamental held as in “quintessence” 
theories [§. 

In particle physics, the first unified theory of the weak and electromagnetic 
interactions @) had as its fundamental bosons just the carriers of the charged 
weak interactions W + , W~, and the photon 7. The next such theory H) had 
a slightly more complicated pattern of gauge bosons — a triplet plus a singlet, 
out of which came not only W + , W~, and 7, but also the neutral weak boson 
Z° , and correspondingly an extra free parameter, the “Weinberg angle.” It 
was of course this latter SU(2)xU(l) theory which has now become part of the 
Standard Model of particle physics. During the early 1970s, however, when the 
experimental data were just becoming available and some of the data appeared 
to contradict the SU(2)xU(l) theory, many other more complicated theories 
were proposed, even by Weinberg 0 , but all these more complicated theories 
ultimately fell by the wayside. 


The development of theories of dark matter may follow a similar pattern. 
By the late 1970s it was becoming clear both that a great deal of dark mat¬ 
ter exists and that the cosmic microwave background (CMB) fluctuation 
amplitude is smaller than that predicted in a baryonic universe. The first non- 
baryonic dark matter candidate to be investigated in detail was light neutrinos 
what we now call “hot dark matter” (HDM). This dark matter is called 
“hot” because at one year after the big bang, when the horizon first encom¬ 
passed the amount of matter in a large galaxy like our own (about 10 12 Mq) 
and the temperature was about 1 keV l§P, neutrinos with masses in the eV 
range would have been highly relativistic. 

It is hardly surprising that HDM was worked out first. Neutrinos were 
known to exist, after all, and an experiment in Moscow that had measured a 
mass for the electron neutrino m(v e ) « 20 eV (corresponding to fi m ~ 1 if 
h were as small as ~ 0.5, since = m(v e )(92h 2 eV)~ 1 ) had motivated espe¬ 
cially Zel’dovich and his colleagues to work out the implications of HDM with 
a Zel’dovich spectrum ( P p (k ) = Ak n with n = 1) of adiabatic primordial fluc¬ 
tuations. But improved experiments subsequently have only produced upper 
limits for m{v e ), currently about 3 eV 0, and the predictions of the adiabatic 
HDM model are clearly inconsistent with the observed universe 0 0 ). 

Cold Dark Matter (CDM) was worked out as the problems with HDM 
were beginning to become clear. CDM assumes that the dark matter is mostly 
cold — i.e., with negligible thermal velocities in the early universe, either 
because the dark matter particles are weakly interacting massive particles 
(WIMPs) with mass ~ 10 2 GeV, or alternatively because they are produced 
without a thermal distribution of velocities, as is the case with axions. The 
CDM theory also assumes, like HDM, that the fluctuations in the dark matter 
have a nearly Zel’dovich spectrum of adiabatic fluctuations. Considering that 
the CDM model of structure formation in the universe was proposed almost 
twenty years ago 0 0’ 0), its successes are nothing short of amazing. As 
I will discuss, the ACDM variant of CDM with fl m = 1 — Ha ps 0.3 appears 
to be in good agreement with the available data on large scales. Issues that 
have arisen on smaller scales, such as the centers of dark matter halos and the 
numbers of small satellites, have prompted people to propose a wide variety 
of alternatives to CDM, such as self-interacting dark matter (SIDM) 0). It 
remains to be seen whether such alternative theories with extra parameters ac- 


tually turn out to be in better agreement with data. As I will discuss below, it 
now appears that SIDM is probably ruled out, while the small-scale predictions 
of CDM may be in better agreement with the latest data than appeared to be 
the case as recently as a year ago. 

In the next section I will briefly review the current observations and the 
successes of ACDM on large scales, and then I will discuss the possible problems 
on small scales. 


2 Cosmological Parameters and Observations on Large Scales 


The table belowQ summarizes the current observational information about the 
cosmological parameters, with estimated la errors. The quantities in brackets 
have been deduced using at least some of the ACDM assumptions. Is is appar¬ 
ent that there is impressive agreement between the values of the parameters 
determined by various methods, including those based on ACDM. In particular, 
(A) several different approaches (some of which are discussed further below) all 
suggest that « 0.3; (B) the location of the first acoustic peak in the CMB 
angular anisotropy power spectrum, now very well determined independently 
by the BOOMERANG 0) and MAXIMAl I§) balloon data HI’ @) and by 
the DASI interferometer at the South Pole HI), implies that Q m + Oa ~ 1; 
and (C) the data on supernovae of Type la (SNIa) at redshifts 2 = 0.4 — 1.2 
from two independent groups imply that Ha — s=s i. Any two of these three 
results then imply that Oa ~ 0.7. The lcr errors in these determinations are 
about 0.1. 

Questions have been raised about the reliability of the high-redshift SNIa 
results, especially the possibilities that the SNIa properties at high redshift 
might not be sufficiently similar to those nearby to use them as standard can¬ 
dles, and that there might be “grey” dust (which would make the SNIa dimmer 
but not change their colors). Although the available evidence disfavors these 
possibilities,^] additional observations are needed on SNIa at high redshift, both 

2 Further discussion and references are given in HU). 

3 For example, SNIa at z = 1.2 and ~ 1.7 apparently have the brightness 
expected in a ACDM cosmology but are brighter than would be expected with 
grey dust, and the infrared brightness of a nearer SNIa is also inconsistent with 
grey dust 
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Table 1: Cosmological Parameters [results assuming ACDM in brackets] 


H 0 = 100 h km s" 1 Mpc” 1 , h = 0.7 ± 0.08 
to = 13 ± 2 Gyr (from globular clusters) 

= [12 ± 2 Gyr from expansion age, ACDM model] 
fl b = (0.039 ±0.006)/if 0 2 (from D/H) 

> [0.035/if o 2 from Lya forest opacity] 

D m = 0.4 ± 0.2 (from cluster baryons) 

= [0.34 ±0.1 from Lya forest P(k)] 

= [0.4 ± 0.2 from cluster evolution] 

~ |Da — 1 ± 1 from SN la 
D m ± Da = 1.02 ± 0.06 (from CMB peak location) 

Da = 0.73 ± 0.08 (from previous two lines) 

< 0.73 (2er) from radio QSO lensing 
D„ £ 0.001 (from SuperKamiokande data) 

;£ [0.1 in ACDM-type models] 


to control systematic effects and to see whether the dark energy is just a cos¬ 
mological constant or is perhaps instead changing with redshift as expected 
in “quintessence” models 0. Such data could be obtained by the proposed 
SuperNova Acceleration Probe (SNAP) satellite H), whose Gigapixel opti¬ 
cal camera and other instruments would also produce much other useful data. 
But it is important to appreciate that, independently of (C) SNIa, (A) cluster 
and other evidence for D m ss 0.3, together with (B) ~ 1° CMB evidence for 
D m ± Da ss 1, imply that Da « 0.7. 

All methods for determining the Hubble parameter now give compati¬ 
ble results, confirming our confidence that this crucial parameter has now been 
measured robustly to a ler accuracy of about 10%. The final result from the 
Hubble Key Project on the Extragalactic Distance Scale is 72±8kms _1 Mpc -1 , 
or h = 0.72 ± 0.08, where the stated error is dominated by one systematic 
uncertainty, the distance to the Large Magellanic Cloud (used to calibrate the 
Cepheid period-luminosity relationship). The most accurate of the direct meth¬ 
ods for measuring distances d to distant objects, giving the Hubble parameter 
directly as Hq = d/v where the velocity is determined by the redshift, are (1) 
time delays between luminosity variations in different gravitationally lensed 
images of distant quasars, giving h ~ 0.65, and (2) the Sunyaev-Zel’dovich 




effect (Compton scattering of the CMB by the hot electrons in clusters of 
galaxies), giving h « 0.63 iH El). For the rest of this article, I will take 
h = 0.7 whenever I need to use an explicit value, and express results in terms 
of ft .70 = Ho/70kms -1 Mpc -1 . 

For a ACDM universe with Q, m = (0.2)0.3(0.4,0.5), the expansion age is 
t 0 = (15.0)13.47(12.41,11.61 )hf 0 1 Gyr. Thus for Q m « 0.3 - 0.4 and h « 0.7, 
there is excellent agreement with the latest estimates of the ages of the oldest 
globular cluster stars in the Milky Way, both from their Main Sequence turnoff 
luminosities El) , giving 12 — 13 ± 2 Gyr, and using the thorium and uranium 
radioactive decay chronometers El), giving 14 ± 3 Gyr and 12.5 ± 3 Gyr, 
respectively. 

The simplest and clearest argument that fl m ss 1/3 comes from compar¬ 
ing the baryon abundance in clusters fb = Mh/M*^ to that in the universe as 
a whole as emphasized by White et al. Since clusters are evi¬ 

dently formed from the gravitational collapse of a region of radius ~ 10 Mpc, 
they should represent a fair sample of both baryons and dark matter. This is 
confirmed in CDM simulations The fair sample hypothesis implies that 


n 


m 





(i) 


We can use this to determine S2 m using the baryon abundance flbh 2 = 0.019 ± 
0.0024 (95% C.L.) from the measurement of the deuterium abundance in high- 


redshift Lyman limit systems 


Using X-ray data from an X-ray flux 


limited sample of clusters to estimate the baryon fraction fb = 0.075 h~ 3 ^ 2 
gives El) Q m = 0.25h -1 / 2 = 0.3 ± 0.1 (using h = 0.70 ± 0.08). Estimating 
the baryon fraction using Suiwaev-Zel’dovich measurements of a sample of 18 
clusters gives fb = 0.077 h~ l , and implies f2 m = 0.25h -1 = 0.36 ± 0.1. 

There is another way to use clusters to measure f2 m , which takes advan¬ 
tage of the fact that the redshift at which structures form depends strongly 
on fl m . This happens because in a low-density universe the growth rate of 
fluctuations slows when, on the right hand side of the Friedmann equation, 


h2 = 8irGp 


k_ 

B 2 


A 
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( 2 ) 


the first (matter) term becomes smaller than either the second (curvature) term 
(for the case of an open universe) or the third (cosmological constant) term. 






As I have already discussed, the A term appears to be dominant now; note that 
if we evaluate the Friedmann equation at the present epoch and divide both 
sides by Hq , the resulting equation is just 

1 = kl m + flfc + Ha • (3) 

Therefore, if we normalize the fluctuation power spectrum P(k) for an f l m = 1 
(Einstein-de Sitter) cosmology and for a ACDM one by choosing erg so that each 
is consistent with COBE and has the same abundance of clusters today, then 
at higher redshifts the low-H m universe will have a higher comoving number 
density of clusters. Probably the most reliable way of comparing clusters nearby 
with those at higher redshift uses the cluster X-ray temperatures; the latest 
results, comparing 14 clusters at an average redshift of 0.38 with 25 nearby 
clusters, give f l m = 0.44±0.12 HI HI). There is greater leverage in this test if 
one can use higher redshift clusters, but the challenge is to find large samples 
with well understood cluster selection and properties. The largest such sample 
now available is from the Las Companas Distant Cluster Survey, which goes 
to redshifts ~ 1, from which the preliminary result is il rn = 0.30 ± 0.12 (90% 

CL) ED. 

3 Further Successes of ACDM 

We have already seen that ACDM correctly predicts the abundances of clus¬ 
ters nearby and at zs, 1 within the current uncertainties in the values of the 
parameters. It is even consistent with P(k) from the Lya forest HI) and from 
CMB anisotropies. Low-fl m CDM predicts that the amplitude of the power 
spectrum P(k) is rather large for k s0.02h/Mpc _ , i.e. on size scales larger 
(k smaller) than the peak in P(k). The largest-scale surveys, 2dF and SDSS, 
should be able to measure P(k) on these scales and test this crucial prediction 
soon; preliminary results are encouraging H^) . 

The hierarchical structure formation which is inherent in CDM already 
explains why most stars are in big galaxies like the Milky Way HU) : smaller 
galaxies merge to form these larger ones, but the gas in still larger structures 
takes too long to cool to form still larger galaxies, so these larger structures 
— the largest bound systems in the universe — become groups and clusters 
instead of galaxies. 


What about the more detailed predictions of ACDM, for example on the 
spatial distribution of galaxies. On large scales, there appears to be a pretty 
good match. In order to investigate such questions quantitatively on the smaller 
scales where the best data is available it is essential to do N-body simulations, 
since the mass fluctuations Sp/p are nonlinear on the few-Mpc scales that are 
relevant. My colleagues and I were initially concerned that ACDM would fail 
this test, since the dark matter power spectrum Pdm(k) in ACDM, and its 
Fourier transform the correlation function ^ m (r), are seriously in disagreement 
with the galaxy data P g {k) and £ g (r). One way of describing this is to say that 
scale-dependent antibiasing is required for ACDM to agree with observations. 
That is, the bias parameter b(r) = [Cg( r )/?dm(?')] 1 ^ 2 , which is about unit y o n 
large scales, must decrease to less than 1/2 on scales of a few Mpc IH’ 

This was the opposite of what was expected: galaxies were generally thought 
to be more correlated than the dark matter on small scales. However, when it 
became possible to do simulations of sufficiently hi gh r esolution to identify the 
dark matter halos that would host visible galaxies M m , it turned out that 
their correlation function is essentially identical with that of observed galaxies! 
This is illustrated in Fig. |l]. 

Jim Peebles, who largely initiated the study of galaxy correlations and 
first showed that £ g (r) ~ (r/ro) -1 ' 8 with ro « 5/i _1 Mpc thought that this 
simple power law must be telling us something fundamental about cosmology. 
However, it now appears that the power law £ g arises because of a coincidence - 
an interplay between the non-power-law (see Fig. [I]) and the decreasing 

survival probability of dark matter halos in dense regions because of their 
destruction and merging. But the essential lesson is that ACDM correctly 
predicts the observed £ g (r). 

The same theory also predicts the number density of galaxies. Using the 
observed correlations between galaxy luminosity and internal velocity, known 
as the Tully-Fisher and Faber-Jackson relations for spiral and elliptical galaxies 
respectively, it is possible to convert observed galaxy luminosity functions into 
approximate galaxy velocity functions, which describe the number of galaxies 
per unit volume as a function of their internal velocity. The velocity function 
of dark matter halos is robustly predicted by N-body simulations for CDM- 
type theories, but to connect it with the observed internal velocities of bright 
galaxies it is necessary to correct for the infall of the baryons in these galaxies 



45) 


which must have happened to create their bright centers and disks. 
When we did this it appeared that ACDM with f l m = 0.3 predicts perhaps too 
many dark halos compared with the number of observed galaxies with internal 
rotation velocities V « 200km s -1 El’ EH). While the latest results from the 



Figure 1: Bottom panel: Comparison of the halo correlation function in an 
ACDM simulation with the correlation function of the APM galaxies HI). The 
dotted curve shows the dark matter correlation function. Results for halos with 
maximum circular velocity larger than 120km s - , 150km s _ , and 200km s -1 
are presented by the solid , dot-dashed , and dashed curves, respectively. Note 
that at scales ^0.3/i _1 Mpc the halo correlation function does not depend on 
the limit in the maximum circular velocity. Top panel: Dependence of bias on 
scale and maximum circular velocity. The curve labeling is the same as in the 
bottom panel, except that the dotted curve now represents the bias of halos 
with V m ax > lOOkms -1 . From Colin et al. EH). 







big surveys now underway appear to be in better agreement with these ACDM 
predictions El’ EH), this is an important issue that is being investigated in 
detail El). 

The problem just mentioned of accounting for baryonic infall is just one 
example of the hydrodynamical phenomena that must be taken into account 
in order to make realistic predictions of galaxy properties in cosmological the¬ 
ories. Unfortunately, the crucial processes of especially star formation and 
supernova feedback are not yet well enough understood to allow reliable cal¬ 
culations. Therefore, rather than trying to understand galaxy formation from 
full-scale hydrodynamic simulations (for example EH) ), more progress has been 
made via the simpler approach of semi-analytic modelling of galaxy formation 
(initiated by White and Frenk EH EH’ EH), recently reviewed and extended 
by Rachel Somerville and me EH)). The computational efficiency of SAMs 
permits detailed exploration of the effects of the cosmological parameters, as 
well as the parameters that control star formation and supernova feedback. We 
have shown EH) that both flat and open CDM-type models with = 0.3 —0.5 
predict galaxy luminosity functions and Tully-Fisher relations that are in good 
agreement with observations. Including the effects of (proto-)galaxy interac¬ 
tions at high redshift in SAMs allows us to account for the observed properties 
of high-redshift galaxies, but only for fl m « 0.3 — 0.5 EH) . Models with Cl m = 1 
and realistic power spectra produce far too few galaxies at high redshift, essen¬ 
tially because of the fluctuation growth rate argument mentioned above. 

In order to tell whether ACDM accounts in detail for galaxy properties, 
it is essential to model the dark halos accurately. The Navarro-Frenk-White 
(NFW) EH) density profile pnfw ( t ) °c r -1 (r+r s ) -2 is a good representation of 
typical dark matter halos of galactic mass, except possibly in their very centers 
(§4). Comparing simulations of the same halo with numbers of particles ranging 
from ~ 10 3 to ~ 10 6 , my colleagues and I have also shown EH) that r s , the 
radius where the log-slope is -2, can be determined accurately for halos with 
as few as ~ 10 3 particles. Based on a study of thousands of halos at many 
redshifts in an Adaptive Refinement Tree (ART) EH) simulation of the ACDM 
cosmology, we tUP found that the concentration c v i r = R v ir/ls_ has a log¬ 


normal distribution, with ler A(logc v i r ) = 0.14 at a given mass 


62[ |63|) 


This 


scatter in concentration results in a scatter in maximum rotation velocities 
of AVjnax/Vmax = 0.12; thus the distribution of halo concentrations has as 




large an effect on galaxy rotation curves shapes as the well-known log-normal 
distribution of halo spin parameters A. Frank van den Bosch HI showed, based 
on a semi-analytic model for galaxy formation including the NFW profile and 
supernova feedback, that the spread in A mainly results in movement along 
the Tully-Fisher line, while the spread in concentration results in dispersion 
perpendicular to the Tully-Fisher relation. Remarkably, he found that the 
dispersion in ACDM halo concentrations produces a Tully-Fisher dispersion 
that is consistent with the observed one.j] 

4 Halo Centers 

Already in the early 1990s, high resolution simulations of individual galaxy 
halos in CDM were finding p(r) ~ r~ a with a ~ 1. This behavior implies that 
the rotation velocity at the centers of galaxies should increase as r 1 / 2 , but the 
data, especially that on dark-matter-dominated dwarf galaxies, instead showed 
a linear increase with radius, corresponding to roughly constant density in the 
centers of galaxies. This disagreement of theory with data led to concern that 
CDM might be in serious trouble EH’ EH) . 

Subsequently, NFW EH) found that halos in all variants of CDM are well 
fit by the pNFw(r) given above, while Moore’s group proposed an alternative 
Pm{?) oc r _3 / 2 (r + r M )~ 3/2 based on a small number of very-high-resolution 
simulations of individual halos EH’ EH> HI . Klypin and collaborators (includ¬ 
ing me) initially claimed that typical CDM halos have shallow inner profiles 
with a ~ 0.2 E3), but we subsequently realized that the convergence tests 
that we had performed on these simulations were inadequate. We now have 
simulated a small number of galaxy-size halos with very high resolution H), 
and find that they range between pnfw and Pm- Actually, these two analytic 
density profiles are almost indistinguishable unless galaxies are probed at scales 

4 Actually, this was the case with the dispersion in concentration 
A (log c v i r ) = 0.1 found for relaxed halos by Jing 1HI, while we HI found 
the larger dispersion mentioned above. However Risa Wechsler, in her disser¬ 
tation research with me HI , found that the dispersion in the concentration at 
fixed mass of the halos that have not had a major merger since redshift z = 2 
(and could thus host a spiral galaxy) is consistent with that found by Jing. We 
also found that the median and dispersion of halo concentration as a function of 
mass and redshift are explained by the spread in halo mass accretion histories. 



smaller than about 1 kpc, which is difficult but sometimes possible. 

Meanwhile, the observational situation is improving. The rotation curves 
of dark matter dominated low surface brightness (LSB) galaxies were measured 
with radio telescopes during the 1990s, and the rotation velocity was typically 
found to rise linearly at their centers Z2 O- EH. But a group led by van 
den Bosch III) showed that in many cases the large beam size of the radio 
telescopes did not adequately resolve the inner parts of the rotation curves, 
and they concluded that after correcting for beam smearing the data are on 
the whole consistent with expectations from CDM. Similar conclusions were 
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Figure 2: High resolution Ha rotation curves (filled circles, solid lines) and HI 
rotation curves for the same galaxies (open circles, dotted lines) from Ref. 64. 
The horizontal bar shows the FWHM beam size of the HI observations. From 
Swaters, Madore, and Trewhella IZil) . 
















reached for dwarf galaxies 1 L3 . Swaters and collaborators showed that optical 
(Ha) rotation curves of some of the LSB galaxies rose significantly faster than 
the radio (HI) data on these same galaxies EH) (see Fig. 2), and these rotation 
curves (except for F568-3) appear to be more consistent with NFW EH). At 
a conference in March 2000 at the Institute for Theoretical Physics in Santa 
Barbara, Swaters also showed a Ha rotation curve for the nearby dwarf galaxy 
DD0154, which had long been considered to be a problem for CDM EH EH); 
but the new, higher-resolution data appeared consistent with an inner density 
profile a « l|] 

Very recently, a large set of high-resolution optical rotation curves has 
been analyzed for LSB galaxies, including many new observations ElD. The 
first conclusion that I reach in looking at the density profiles presented is that 
the NFW profile often appears to be a good fit down to about 1 kpc. However, 
some of these galaxies appear to have shallower density profiles at smaller radii. 
Of the 48 cases presented (representing 47 galaxies, since two different data sets 
are shown for F568-3), in a quarter of the cases the data do not probe inside 1 
kpc, and in many of the remaining cases the resolution is not really adequate 
for definite conclusions, or the interpretation is complicated by the fact that 
the galaxies are nearly edge-on. Of the dozen cases where the inner profile is 
adequately probed, about half appear to be roughly consistent with the cuspy 
NFW profile (with fit a^0.5), while half are shallower. This is not necessarily 
inconsistent with CDM, since observational biases such as seeing and slight 
misalignment of the slit lead to shallower profiles EH) . Perhaps it is significant 
that the cases where the innermost data points have the smallest errors are 
cuspier. 

I think that this data set may be consistent with an inner density pro¬ 
file a ~ 1 but probably not steeper, so it is definitely inconsistent with the 
claims of the Moore group that a M.5. But very recent work by Navarro and 
collaborators EH) has shown that Moore’s simulations did not have adequate 
resolution to support their claimed steep central cusp; the highest-resolution 
simulations appear to be consistent with NFW, or even shallower with a « 0.75. 
Further simulations and observations, including measurement of CO rotation 
curves EH), may help to clarify the nature of the dark matter. 


5 Swaters (private communication) and Hoffman have subsequently con¬ 
firmed this with better data, which they are preparing for publication. 



It is something of a scandal that, after all these years of simulating dark 
matter halos, we still do not have a quantative — or even a qualitative — theory 
explaining their radial density profiles. In her dissertation research HI), Risa 
Wechsler found that the central density profile and the value of r s are typically 
established during the early, rapidly merging phase of halo evolution, and that, 
during the usually slower mass accretion afterward, r s changes little. The 
mass added on the halo periphery increases R v i r , and thus the concentration 
c v i r = R v i r /r s . Now we want to understand this analytically. Earlier attempts 
to model the result of sequences of mergers (e.g., EH’ led to density 
profiles that depend strongly on the power spectrum of initial fluctuations, in 
conflict with simulations (e.g. EH)). Perhaps it will be possible to improve on 
the simple analytic model o f ma ss loss due to tidal stripping during satellite 
inspiral that we presented in ). Avishai Dekel and his students have recently 
shown that including the tidal puffing up of the inspiralling satellite before tidal 
stripping can perhaps account for the origin of the cusp seen in dissipationless 
simulations, independent of the power spectrum. They argue that the profile 
must be steeper than a = 1 as long as enough satellites make it into the halo 
inner regions, simply because for flatter profiles the tidal force causes dilation 
rather than stripping. The proper modeling of the puffing and stripping in 
the merger process of CDM halos may also provide a theoretical framework 
for understanding the observed flat cores as a result of gas processes; work on 
this by Ari Mailer and Dekel is in progress. Reionization and feedback into the 
baryonic component of small satellites would make their cores puff up before 
merging. This could cause them to be torn apart before they penetrate into 
the halo centers, and thus allow a < 1 cores. 

Another possible explanation for flatter central density profiles involving 
the baryonic component in galaxies has recently been proposed 0) , in which 
the baryons form a bar that transfers angular momentum into the inner parts 
of the halo. It is not clear, however, that this effect could be very important in 
dark matter dominated dwarf and LSB galaxies that have small or nonexistent 
bulge components. 

It would be interesting to see whether CDM can give a consistent account 
of the distribution of matter near the centers of big galaxies, but this is not 
easy to test. One might think that big bright galaxies like the Milky Way 
could help to test the predicted CDM profile, but the centers of such galaxies 




are dominated by ordinary matter (stars) rather than dark matter.^ 

5 Too Much Substructure? 


Another concern is that there are more dark halos in CDM simulations with 
circular velocity 14^30 km s -1 than there are low- V c galaxies in the Local 
Group H' 0. A natural solution to this problem was proposed by Bullock 
et al. 0), who pointed out that gas will not be able to cool in V c s30 km s -1 
dark matter halos that collapse after the epoch of reionization, which occured 
perhaps at redshift z re ion ~ 6 0). When this is taken into account, the 
predicted number of small satellite galaxies in the Local Group is in good 
agreement with observations 0 0). It is important to develop and test this 
idea further, and this is being done by James Bullock and by Rachel Somerville 
and their collaborators: the results to date (e.g. 0 0)) look rather promising. 

"0 ) now agree that astrophysical effects will keep 
As a result, theories such as warm dark matter 


Other groups (e.g. 
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most of the subhalos dark. 

(WDM), which solve the supposed problem of too many satellites by decreasing 
the amount of small scale power, may end up predicting too few satellites when 
reionization and other astrophysical effects are taken into account 0) . 

The fact that high-resolution CDM simulations of galaxy-mass halos are 
full of subhalos has also led to concerns that all this substructure could prevent 
the resulting astrophysical objects from looking like actual galaxies 0). In 
particular, it is known that interaction with massive satellites can thicken or 
damage the thin stellar disks that are characteristic of spiral galaxies, after 
the disks have formed by dissipative gas processes. Ho wever, detailed simula¬ 


tions 


have shown that simpler calculations Ii25|) had overestimated the 
extent to which small satellites could damage galactic disks. Only interaction 
with large satellites like the Large Magellanic Cloud could do serious damage. 
But the number of LMC-size and larger satellites is in good agreement with the 
number of predicted halos 0, which suggests that preventing disk damage 
will not lead to a separate constraint on halo substructure. 


6 Navarro and Steinmetz had claimed that the Milky Way is inconsistent 
with the NFW profile 0) , but they have now shown that ACDM simulations 
with a proper fluctuation spectrum are actually consistent with the data 0). 




6 Angular Momentum Problems 


As part of James Bullock’s dissertation research, we found that the distribu¬ 
tion of specific angular momentum in dark matter halos has a universal pro¬ 
file ). But if the baryons have the same angular momentum distribution as 
the dark matter, this implies that there is too much baryonic material with low 
angul ar momen tum to form the observed rotationaUv supported exponential 
disks 101 > 102). t^g | on g i 3 een assumed (e.g. El’ 01 ) that the baryons and 
dark matter in a halo start with a similar distribution, based on the idea that 
angular momentum arising from large-scale tidal torques will be similar across 
the entire halo. But as my colleagues and I argued recently, a ke y implication of 
our new picture of angular momentum growth by merging 110) j s that the DM 
and baryons will get different angular momentum distributions. For example, 
the lower density gas will be stripped by pressure and tidal forces from infalling 
satellites, and in big mergers the gaseous disks will partly become tidal tails. 
Feedback is also likely to play an important role, and Mailer and Dekel (in 
preparation) have shown using a simple model that this can account for dat a 
on the angular momentum distribution in low surface brightness galaxies 
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A related concern is that high-resolution hydrodynamical simulations of 
galaxy formation lead to disks that are much too small, evidently because 
formation of baryonic substructure leads to too much transf er of angular mo¬ 
mentum and energy from the baryons to the dark matter ft^) But if gas 
cooling is inhibited in the early unive rse, m ore realistic disks form EEH) , more 


so in ACDM than in = 1 CDM 


107) 


Hydrodynamical simulations also 


indicate that this disk angular momentum problem may be resolved if small 
scal e power is suppressed because the dark matter is warm rather than cold 


12 ^), which I discuss next. 


7 Alternatives to ACDM? 

Because of the concerns just mentioned that CDM may predict higher densities 
and more substructure on small scales than is observed, many people have 
proposed alternatives to CDM. Two of these ideas that have been studied in 
the greatest detail are self-interacting dark matter (SIDM) EH) and warm dark 
matter (WDM). 

Cold dark matter assumes that the dark matter particles have only weak 










interactions with each other and with other particles. SIDM assumes that the 
dark matter particles have strong elastic scattering cross sections, but negligi¬ 
ble annihilation or dissipation. The hope was that SIDM might suppress the 
formation of the dense central regions of dark matter halos, although the large 
cross sections might also lead to high thermal cond uctiv ity which drains energy 
from halo centers and could lead to core collapse 222, and which also causes 
evaporation of galaxy halos in clu ster s, resulting in violation of the observed 
“fundamental plane” correlations 2L2) . But in any case, self-interaction cross 


sections large enough to have a significant effect on the centers of ga l axy- mass 


halos will make the ce nters of g alaxy clusters more spherical ill’ LL2 ) and 


perhaps also less dense 112’ 111) than gravitational lensing observations 111) 
indicate. 

Warm dark matter arises in particle physics theories in which the dark 
matter particles have rela tively high thermal velocities, for example because 
their mass is si keV Iii2l), comparable to the temperature about a year after 
the Big Bang when the horizon first encompassed the amount of dark matter 
in a large galaxy. Such a velocity distribution can suppress the formation of 
structure on small scales. Indeed, this leads to constraints on how low the 
WDM particle mass can be. From the requirement that there is enough small- 
scale power in the linear power spectrum to reproduce the observed properties 
of the Lya forest in quasar spectra, it follows that this mass must exceed 
about 0.75 keV 


117) 


The requirement that there be enough small halos 
to host early galaxies to produce the floor in metallicity observed in the Lya 
forest systems, and early galaxies and quasars to reionize the universe, pro babl y 
implies a st rong e r low er limit on the WDM mass of at least 1 keV 
Simulations 


118) 


120 ) 


do show that there will be far fewer small satellite halos 
with AWDM than ACDM. However, as I have already mentioned, inclusion of 
the effects of reionization may make the observed numbers of satellite galaxies 
consistent with the predictions of ACDM EB), in which case AWDM may 
predict t oo few sm all satellite galaxies. Lensing can be used to look for these 
subhalos 


121, pa) 


and m ay already indicate that there are more of them than 
expected in AWDM EH ). Thus it appears likely that WDM does not solve all 
the problems it was invoked to solve, and may create new problems. Moreover, 
even with an initial power spectrum truncated on small scales, simulations 
appear to indicate that dark matter halos nevertheless have density profiles 
















much like those in CDM 
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I 


(although doubts have been ex press ed 


about the reliability of such simulations because of numerical relaxation ll^l)). 


But WDM does lead to lower co ncen t ratio n halos in better agreement with 
observed rotation velocity curves 


126, 127) 


One theoretical direction that does appear very much worth investigating 
is A CDM with a tilt n ~ 0.9 in the primordial power spectrum P p (k ) oc k n 


. Such tACDM cosmology is favored by recent measurements of the power 


spectrum of the Ly a forest £§) and appears to b e con sistent with the latest 
CMB measurements and all other available data Our simple analytic 

model 0) predicts that the concentration of halos in tACDM will be approx¬ 
imately half that in LCDM, which appears to be true in a trial simulation by 
A. Kravtsov. While this does not resolve the cusp problem, it is a step in the 
right direction which may lessen the conflict with galaxy rotation curves. 


8 Outlook 

The successes of the CDM paradigm are remarkable. Except possibly for the 
density profiles at the centers of dwarf and low surface brightness galaxies, the 
predictions of ACDM appear to be in good agreement with the available ob¬ 
servations. The disagreements between predictions and data at galaxy centers 
appear to occur on smaller scales than was once thought, but as the data im¬ 
prove it is possible that the discrepancies on Si kpc scales may ultimately show 
that CDM cannot be the correct theory of structure formation. However, it 
appears to be better than any alternative theory that has so far been studied, 
even though these alternative theories have more adjustable parameters. 

This article started by discussing the analogy between the effort to un¬ 
derstand dark matter and structure formation in modern cosmology and the 
effort to understand particle physics in the 1960s and 1970s. In both cases, 
the result was a “standard model” which has guided further work and led to 
great progress in both theory and observation/experiment. But in both cases, 
the standard model is not an ultimate theory, and the search is on for a better 
theory. In the case of particle physics, there is a leading candidate: super- 
symmetry, and perhaps ultimately string or M theory. Here the analogy fails, 
because I am not aware of any theory that has all the virtues of CDM but 
which avoids its possible failure at the centers of galaxies. The quest for such 
a theory is a worthwhile goal. But for many purposes, including studies of the 











formation and evolution of galaxies and their large scale distribution, the CDM 
standard model may still remain very useful. And maybe it is even true. 
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